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Abstract 

We propose a possible mechanism of pseudo spin gap anomaly(PSGA) in 
magnetically coupled bi-layer cuprates without any fermion pairing instability. 
In our proposal PSGA does not necessarily require the spin-charge separation or 
the breakdown of the Fermi liquid description of a normal state of the cuprate 
superconductor. 

The low energy magnetic excitations are mainly governed by the itinerant 
nature of the intra-layer system and the inter-layer antiferromagnetic coupling. 
No matter how weak the bare inter-layer coupling is, it can be dramatically en- 
hanced due to the intra-layer spin fluctuations. As the temperature decreases 
near the antiferromagnetic phase boundary the strongly enhanced inter-layer 
correlation induces the inter-layer particle-hole exchange scattering processes 
that tend to enhance the inter-layer spin singlet formation and kill the triplet 
formation. We propose that the coupling of spin fluctuations on the adjacend 
layers via the strong repulsive interaction between parallel spins travelling on 
each layer give rise to the dynamical screening effects. As a result the low energy 
part of the spin excitation spectrum is strongly suppressed as the temperature 
decreases near the antiferromagnetic phase boundary. We ascribe PSGA to this 
dynamical screening effects. 
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1 Introduction 



Since the discovery of the high Tc cuprate superconductor in 1986, puzzhng features 
in their normal state due to strong correlation especially the spin dynamics of 
them Pi have provoked a great deal of controversy. In the high superconductor, a 
superconducting phase lies just near the Mott insulating phase where only the spin 
degree of freedom survives. Thus even in a metallic phase the antiferromagnetic spin 
fluctuations play essential roles to understand the physical nature of these materials. 
Especially the spin dynamics of bi-layer cuprates have recently attracted a lot of 
attention, since the so-called pseudo spin gap anomaly(PSGA), that is one of the 
most serious problem in this field, seems to be peculiar to bi-layer cuprates. 

PSGA was first observed by Yasuoka^ as the phenomena that NMR longitudinal 
relaxation rates, 1/TiT deviate from Curie- Weiss like temperature dependence well 
above the superconducting transition temperature and show a broad peak around the 
so-called spin gap temperature before it decreases smoothly down through Tc with 
slight change of slope. Up to now the clear signature of PSGA in 1/TiT can be 
observed in YBa2Cu306+.. § |, B |, 0, B i and , YBa2Cu408|T2|, |13|, 0, 
Y2Ba4Cu70i5[jl5|, LaBa2Cu30y|T^ , and Bi2Sr2CaCu20t,|lT7|, |18| that are all 'bi-layer' 



cuprates where there are two equivalent Cu02 layer in a unit cell. Furthermore in the 
lightly doped YBa2Cu306+a;, the inelastic neutron scattering cross sections deviate 
from a linear u) dependence in the low energy region and are strongly suppressed 
above the superconducting transition temperature [2^ . 

On the other hand in typical 'mono-layer' cuprate La2_xSrj.Cu04, 1/TiT begin 
to deviate slightly from Curie- Weiss law below T ~ 200K||2^, but do not show clear 
peak structure above Tc. As for the neutron experiments on LSCO system, only 
recently an overall energy spectrum of magnetic fluctuations was reported in 
homogeneous single crystal of Lai.85Sro.i5Cu04. They first reported the gap formation 
in LSCO system belowTc and succeeded in detecting the superconducting phase of this 
system. However they didn't detect the pseudo-spin gap anomaly in this compound. 
Furthermore the single layer over-doped cuprate Tl2Ba2CuOj; exhibits a very broad 
peak in both of 1/TiT and I/T2 slightly above Tc. As Yasuoka pointed outfl^, |25|, |26[|, 
this behavior should be atributed to the band effects and cannnot be regarded as the 
spin-gap behavior. These experimental facts indicate that PSGA is peculiar to bi- 
layer cuprates. 

The most important feature of PSGA is that the strong suppression of the low 
energy spin excitations does not accompany the clear gap like structure in the low en- 
ergy charge excitations. By reason of this fact, the spin gap in high Tc superconductor 
is called the pseudo gap. The suppression of the spin excitations clearly indicates that 
the singlet pair formation is favored and as a result the magnetic degree of freedom 
tends to be killed^T^. Therefore to understand PSGA phenomena we have to seek 
some mechanism that favors the singlet pair formation without any suppression of 
the charge excitations. 

From theoretical viewpoints there are two fixed points to describe low energy ex- 
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citations in a normal state of the cuprate superconductor. One fixed point is the 
Fermi hquid fixed point where the low energy excitations are described in terms of 
quasi-paricle excitations. If we try to explain PSGA within this framework, we have 
to seek a mechanism that enhances the singlet pair formation without any pairing in- 
stability since pairing instability in a quasi-particle system means the superconducting 
transition where both of the spin and charge degree of freedom condensate into pair. 
Another fixed point is the non-Fermi liquid picture with spin charge separation where 
the low energy excitations are spin 1/2 chargeless fermions (spinous) and charged spin- 
less bosons(holons)|]l|. In this picture PSGA is ascribed to the condensation of the 
spinous into the singlet pair. 

The Hubbard model is the simplest model to describe the effects of the strong 
correlation within the Fermi liquid picture. In the Hubbard model, we take into 
account the electronic correlations by including the on-site repulsion, U, between 
two electrons with opposite spins. Then the simplest scheme to describe the collec- 
tive spin excitations is the random phase approximation(RPA), namely paramagnon 
theory [Q. Apart from PSGA, quantitative behavior of dynamical spin fluctuations 



can well be described within the simple RPA[p9|, 30, 31]. We can go beyond RPA 



by treating the feedback effects of the spin fluctuations on the irreducible polariza- 



tion in a self consistent manner|32, p3, 34|. More sophisticated treatment based on 



Baym-Kadanoff formalism that treats the exchange of fluctuations in the particle-hole 
and particle-particle channels in a self consistent manner has been proposed 
Phenomenological approach to the antiferromagnetic spin fluctuations proposed by 
Millis, Monien and Pines also lays its microscopic foundation on the Fermi liq- 
uid picture. However it seems difficult to explain PSGA based only on the Hubbard 
model, since the Hubbard interaction never distinguishes the singlet and triplet pair 
selectively. In the Hubbard model there appear only one scattering channel between 
quasi-particles with opposite spins and as a result it cannot describe the tendency of 
the enhanced singlet. 

On the other hand the non-Fermi liquid picture with spin-charge separation stems 
from the t-J model that can be regarded as the strong coupling limit of the Hub- 
bard model ^8|. This model can give magnetic phase diagram within the mean field 
level p9[]. Recently the extention of the t-J model to the bi-layer cuprates was also 



reported]^, ^ In the t-J model the presence of the spin dependent interaction 
is essential. It produces three fundamental scattering channels, one with spin-flip and 
two without spin-flip. These scattering channels can be decomposed into the singlet 
and triplet channel and spin dependent nature of scattering process can well be de- 
scribed. Concerning PSGA there have been some proposals mainly from the gauge 
field theoretical approach based on the t-J model^, |3|, ||, |5|, ||, g^l- Altshuler 
and Ioffe||45|| first pointed out the importance of bi-layer correlation. They claimed 
that if the spins on each layer form a spin liquid state, this transition can be de- 
scribed as a spinon pairing on adjacent layers, that leads to the spin-gap anomaly. 
As for the spinon condensation into the singlet pair, Ubbens and Lee claimed 
that f/(l) gauge fluctuations tend to destroy the intra-layer spinon pairing, although 
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it is less effective in destroying the inter-layer pairing P6[. However recently Wen 



and Lee|^ reported that the spin gap phase is best understood as the intra-layer 
staggered ffux phase where the SU{2) symmetry is preserved away from half filling. 
Therefore within the framework of 2-dimensional gauge theory there still exist some 
ambiguity on whether PSGA is intrinsic to bi-layer structure or not. 

However it remains an unsettled question whether the low energy excitations in 
cuprate can be described as the Fermi liquid or not. In the present paper, we should 
like to explore a different possibility to understand PSGA without any fermion con- 
densation, which has never been examined. As stated above to understand PSGA 
the existence of spin dependent interaction is essentially important. However this 
does not necessary mean the spin-charge separation really takes place. Motivated 
by this consideration, we try to understand PSGA without breakdown of the Fermi 
liquid description. To compensate the drawback of the Hubbard model that it can- 
not distinguish the singlet pair formation selectively, we consider the magnetically 
coupled bi-layer cuprates and try to understand the PSGA without any pairing in- 



stability. Our motivation was first provoked by the idea by loffe et.al\i8\, where the 



strongly developed inter-layer spin fluctuations dramatically enhance the inter-layer 
interaction. In the present paper we treat the intra-layer spin fluctuations based on 
the Hubbard model. As the temperature decreases, the developed intra-layer spin 
fluctuations dramatically enhance the inter-layer magnetic correlation near the mag- 
netic phase boundary. Therefore no matter how weak the bare inter-layer magnetic 
coupling is, it can affect the whole structure of the magnetic excitations of the system 
at the low temperature near the magnetic phase boundary. 

Just near the antiferromagnetic phase boundary, the enhanced inter-layer corre- 
lations strongly induce the inter-layer particle-hole exchange scattering processes. In 
this paper we show that due to these scattering processes the inter-layer singlet forma- 
tion is strongly enhanced while the triplet formation is strongly suppressed. In other 
words interaction between parallel spins becomes strongly repulsive. As a result we 
can ascribe the PSGA to the dynamical screening of the low energy spin excitations 
due to this repulsion. Actually we can show that as the temperature decreases and 
the intra-layer spin fluctuations develop, the imaginary parts of the whole suscepti- 
bility reveal apparent gap-like structure in the low energy region. We can assign the 
PSGA to this mechanism. 

The outline of this paper is as follows. In section 2, we present the model and 
formalism. In section 3, we present how dramatically the inter- layer bare interaction 
is enhanced by the intra-layer spin fluctuations. We pay our attention to the spin 
rotational invariance relation. In section 4, we investigate the inter-layer particle-hole 
exchange scattering processes driven by the enhanced inter-layer intetaction in some 
detail. In section 5, we show that the coupling of spin fluctuations on the adjacend 
layers via the inter-layer processes give rise to the dynamical screening effects. In 
section 6, we present the relations between microscopic spin correlations obtained in 
the present approach and experimentally observed quantities in magnetically coupled 
bi-layer cuprates. In section 7, We present the numerical results on the analytically 
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continued dynamical susceptibilities and show that our proposal can actually be re- 
alized. 

2 Model and Formalism 

In this section we present our model and basic formalism to treat the spin fluctuations 
in magnetically coupled bi-layer. We consider YBa2Cu306+a: (Y123) compounds as 
representatives of the bi-layer family. As shown in Fig. 1. there are a pair of equivalent 
Cu02 layers per unit cell that are separated by the charge reservoir layers. We ignore a 
small orthorhombic distortion of the Cu02 lattice. The distance between the adjacent 
Cu02 layers is approximately c* ~ 3.2rA. On the other hand the distance between 
the pair of layers is approximately c ~ 8.2rA. 

2.1 Model Hamiltonian 

We start with the Hubbard Hamiltonian with small direct inter-layer hopping 

— X/ X/ X/ ^ij^n,j,a^n,i,a + ^-L X/ X/ '^n,i,o-<^m,i,<T + ^ X/ X/ '^n,i,'['nn,i,i: (l) 

where c|j j ^ {cn,i,a) is a creation (annihilation) operator of an itinerant carrier belong- 
ing to the i-th Cu site of the n-th layer (n = 1,2) with the spin projection a and 
^n,i,o- = Cn,i,J Cn,i,a- Furthermore tij denotes the intra-layer hopping integral between 
the i-j sites and denotes the inter-layer hopping integral. As for the intra-layer 
on-site Hubbard repulsion, f/, we assume the intermediate coupling scheme, namely 
U ~ 3i, where t is the hopping integral between the nearest neighbor sites. On the 
other hand, we can safely assume, t± <^ U and therefore we can treat the direct 
inter-layer hopping process in a t-J like manner. Thus we obtain 

n^n// + n±, (2) 

where 

(3) 

i,j an in 

denotes the intra-layer Hubbard Hamiltonian and 

'Hi. = ^-L E E (-'■ ~ ''^n,i,-a)Cn^i,aCm,i,ai^ — nm,i-a) + 2 Jo±X^ X^ Sm,i ' Sn,i, (4) 
a n^m i m^n 

is the t-J like counterpart that gives rise to the inter- layer magnetic coupling. The 
spin fluctuation operator is introduced by 

'^1"'' ^ 9 X! 4i,i,af^ al3Cn,i,p, (5) 
^ ad 

where cr denotes usual Pauli matrices and we set 7;, = 1. We note that the hopping 
term in 7i± depends on the carrier density and near the antiferromagnetic phase 
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boundary it is reduced, while the magnetic couphng is unhkely to be sensitive to the 
doping. Furthermore, as will be shown later, near the phase boundary the inter-layer 
coupling is dramatically enhanced by the intra-layer spin fluctuations. Based on this 
fact, in H_L we keep only the magnetic term. 

Thus we are lead to the assumption that the spin fluctuations in bi-layer cuprates 
can well be described by the effective Hamiltonian 

^// = J2J2J2 Uj4i,j,aCn,i,a + U J^Yl ^n,^,!^",*,^' (6) 
i,j an in 

and 

^m,i ' ^n,i- (7) 

By taking the lattice Fourier transform 



the effective Hamiltonian can be written by 

^// = E E E^n,fccl_fc_,C„_fc_^ 

and 

^± ^ = 2Jo±E E ■ Sni-q). (10) 

where N is the number of the lattice sites and 

denotes the spin fluctuation in the momentum space. Here we note the momentum 
k and q are the two-dimensional momentum since in our model the itinerant carriers 
are confined into a single plane. In the kinetic term, we have included the hopping 
process to the second nearest neighbor sites 

El If. = e2,fc = ~ — 2t(cos k^a + cos kyo) — 2t' cos kr^a cos kyO 

— 2t"(cos 2kxa + cos 2kya), (12) 

where we put t' = —t/5 and t" = t/4 to reproduce the Fermi contour of YBCO 
system [|1|, ||]. 

Now we briefly comment on the experimental evidences for the inter-layer mag- 
netic coupling. It is well established that there remains the strong intra-layer anti- 
ferromagnetic spin-spin coupling which is estimated as Jq// = 80^3q meV from the 
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neutron experiments [^] . On the other hand, although existence of the inter-layer 
antiferromagnetic spin-spin interaction has been widely accepted, question about the 
magnitude of this interaction had been under debates. Concerning this issue, the 
early neutron scattering experiment gave very weak strength of inter-layer magnetic 
coupling Jo± ~ O.O6J0//. This estimation was based on the fact that the expected 
optical gap of spin wave excitation, 2JJq//Jo±, in an insulating phase had never 



been observed even up to 60 meV||51||. However recently it was reproted that mid- 
infrared transmission and reflection mesurements picked up the optical branch at 
178.0 meV [^] and therefore Jo± can be estimated as Jo± ~ 0.55Jo//. Moreover re- 
cent NMR cross-relaxation mesurements of Y2Ba4Cu70i5 |TB] also suggest that the 
inter-layer spin-spin coupling can reach Jq± ~ 0.25Jo//. Tus we can say that the 
importance of the inter-layer magnetic coupling should be duely recognized to under- 
stand the magnetic properties of bi-layer cuprates. 



2.2 Dynamical Spin Susceptibility 

In the case of magnetically coupled bi-layer, the dynamical spin susceptibility can be 
written in the thermal Green's function formalism as 

rl/T 

= / rfre'-'^<T,K(g,r)^f(-q,0)]> 
Jo 

= ^E^^.r(^rn'^(9;^^0<, (13) 

fiuXp 

where 

k,k' 

(14) 

denotes the spin dependent polarization function. Now m,n are layer indices, = 
+, — , z are the spin indices, /i, i/. A, p =t, [ are the spin projection, and ui = 2nTl is 
a bosonic Matsubara frequency. Furthermore 

<,fc,„M=^-P(-^)<,fc,J0)exp(-r7^) (15) 

represent imaginary time dependent creation operator where we set ^ = 1 and kB = 
1. The thermal average < ••• >= Tr(e~^''^ • • •)/Tre~^/"^ is taken under the full 
Hamiltonian. Due to the symmetry under exchange of layers, the spin susceptibility 
Xmn has only two independent components with respect to m, n: the intra-layer 
(diagonal) correlation Xii = X22 = X// ^"^^ the inter-layer [off- diagonal) correlation 
X12 = X21 = xt ■ Furthermore each element has two independent components in 
spin space: transverse counterpart Xmn ^^"^ longitudinal counterpart Xmn' which are 
given by 

(16) 
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Throughout this paper we consider the paramagnetic pahse and therefore the dy- 
namical spin susceptibihties must satisfy the rotational symmetry relation in spin 
space; 

2Xm„(g; i^l) = Xmn(9; ^^l) ■ (17) 

We should always check this relation when we carry the theoretical investigation a 
stage further. Since we assumed there is no inter-layer carrier hopping, the non- 
interacting counterpart can simply be written as 

XotniQ'^ i^i) = 2X0mn(9; i^i) = SmnXoiq; i^Ji), (18) 



where 



T 

i£n k 



(19) 



1 ^fc.g(^)(^fc+g-^fc) 



2N 



k 



Now 



Qo{k;ier 



1 



(20) 



^k 

is a Green's function for a free carrier where = (2n + l)'7rT is a fermionic Matsubara 
frequency. The thermal extinction factor is given by 



^k 



^fc q{T) = tanh(^^) - tanh(^ 



2T 



■2T' 



(21) 



where Cfc = Sfe ~ ^■ 

First we consider the intra-layer spin fluctuations driven by the Hubbard interac- 
tion within the frame work of random phase approximation (RPA). Taking into the 
effects of the interaction, instead of (^), the dynamical susceptibility is rewritten in 
a form, 

T 

X"^{q; = J2 ^o(fc + q, i^n + uji)-f''^{q; iuJi)goik; iSn), (22) 



where "^"^{q) denotes the RPA triangle vertex. Usual RPA for Hubbard model ||28[ 
gives 
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7^^' "{q-.iui) 



1 



1 - Uxo{q]iuJi)' 

1 

1 - U^xo{q;'iu;iy' 

UXoiq;iuJi) 
l-U^X^iq-iujif 



(23) 
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These vertices satisfy the spin-rotational invariance relation 



(24) 



where 'j{q;iui) = [1 — Uxo{q',i^^i)] ""^ is a Stoner factor. These vertices give the 
dynamical spin susceptibilities; 

iuJi) = Xoiq; iuJih^' {q; ^^i)^ 

x'^'i.q] ii^i) = Xoiq; iuJi)Y''iq; ^^i), 
x'^'^'^iq; ii^i) = Xoiq; iuJi)Y'^'' iq; i^i), 



(25) 



r%q;iui) = -J2[x'''''iq;^^i)-x'''-''iq;^^i)]- 

where we used the relation (|T^). We note that in this step the spin rotational sym- 
metry relation 

x^'iq; i^i) = '^x'^iq; i^i) = xiq; i^i), (26) 

is satisfied where we define the RPA spin susceptibility 



xiq; 



xo[q;ii^i) 



1 - Uxoiq;i^i) 



(27) 



3 Enhanced Inter-layer Interaction 



The fundamental processes induced by the inter-layer interaction are given through 
the decomposition 

5i ■ 52 = \iStS2 + S^S^) + SIS^2- (28) 

By rewriting this expression in a second quantized form, we get three scattering 
processes between propagating spins on each layer. The first term of (^) produces 
the scattering process with spin flip (Type-A). The second term of ( pSD produces 
two fundamental processes; scattering processes between parallel spins (Type-B) and 
between anti-parallel spins (Type-C). These processes are shown in Fig. 2. These 
three scattering processes can couple to the intra-layer spin fluctuations in each of 
adjacent layers. In Fig. 3. we show how the inter-layer exchange interaction vertices 
couple to the intra-layer RPA spin fluctuations. Thus the intra-layer spin fluctuations 
on each layer can couple to each other via these scattering vertices as shown in Fig. 4. 
As is shown in Fig. 4 (a), the Type-A scattering can couple only to the intra-layer 
transverse spin fluctuations. Further Type-B and Type-C scattering can couple only 
to the intra-layer longitudinal spin fluctuations, as depicted in Fig. 4(b),(c). Then 
the enhanced inter-layer interaction for each channel can be given respectively by: 



-J| iq;iuji) 



[q;tuJi)—-f^ 



iq;i^i) 



-iiq;i^if 



(29) 
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for the type-A scattering and 

= ^iiq;mr, (30) 

and for the type-B scattering 

- ^r^-^'iqituJi) = -Y^--{q;tuj,yI^Y'--{q;iu;i)+r'--iq-,tcuiy^r^'--iq;iu;i) 
+ r'^iq; iu;i)^r'-''(q; io^i) - l^'^iq; iu;i)^T'''-''(q; m) 



^ 7{q;^<^i)% (31) 

for the type-C scattering. Thus scattering vertices corresponding to each channel are 
enhanced in the same manner. This situation comes from the rotational symmetry 
of the interaction in the spin space. Thus we introduce the enhanced inter-layer 
interaction 

Jj_(g; iuji) = Jo±7(g; i^if- (32) 

We can see directly from these expressions that the inter-layer Heisenberg interaction 
can be strongly enhanced by the strongly enhanced intra-layer spin fluctuations. 



4 Inter-layer Exchange Scattering Processes 

When we consider the spin fluctuating proceeses driven by the strongly enhanced 
inter-layer interaction, we should bear in mind that since in our treatment we neglect 
the inter-layer carrier hopping, the electron-hole pair bubble laid across adjacent 
layers doesn't exist. Therefore the most dominant process driven by the inter-layer 
interction is the electron-hole exchange scattering where an electron and a hole run 
in different layers. This kinds of processes are expected to develop dramatically just 
near the antiferromagnetic phase boundary and modify the simple RPA susceptibility. 
These scattering processes are described in terms of the T-matrices deflned by 

rr'^iK k'; q) = fj^ dre^-^^ < ^.[4^^, ^(r)c,_fe,^^^,(r)c^^^ ;0)c,,fe,,(0)] >, (33) 
where /i, u, A, p =1 or J,. 
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4.1 Scattering T-matrices 

The general expression ( P^D gives spin dependent T-matrices, 

Tr'^\k,k';q)=rl^^^^{k,k';q), 
[ Tl^'^\k,k';q)=Tl^'^^ik,k';q). 



(34) 



To clarify the nature of these spin dependent scattering processes, let us decompose 
them into the singlet and triplet channel. Since in our case T-matrices must satisfy the 
spin rotational symmetry relation, the scattering vertex can be written in a operator 
form in the spin space, 

71 = T^p + Tla^l • (72 (35) 

where cti and (T2 denote the Pauli matrices that represent spin on the layer-1 and 
layer-2 respectively. Then cri ■ (T2 has eigen values 1 and —3 that correspond to the 
triplet and singlet pair respectively. We can decompose T± into the singlet and triplet 
channel as 



triplet 



rrsinglet 



- Tip + TL(T5 

Tip — 371o-- 



(36) 



Thus we get 



r-r-singlet , Qrrtriplet ^triplet 

1 I = : 1 



-triplet 



/j-singlet 



-CTi ■ (72. 



(37) 



4 4 
On the other hand, by noting cri ■ (T2 = 2{afa2 + cricrt) + ^i*^!) '^^ get 

' 'T-Ti.Ti 1 ( q- singlet . q-triplet'^ 

/-7-tTiTT ^triplet 

/T-Ti,J.T 1 (rrtriplet ^singlet\ 

From these equations we can directly check the spin rotational invariance relation 

^n,TT_^_TUT^^n,n_ (39) 



(38) 



TT 11 1 • ^sinqlet i ^zripiez .ji ^fj,u,Ap 

h mally we can combme ij^ and with IJ^ like 



-triplet 



-liv,\p 



rj-singlet 
rj-triplet 



-jTT.TT 



(40) 
(41) 



In Fig. 5 we show the diagrammatic representation for the particel-hole ladder 
processes where the n-th order ladders are shown. As is shown in Fig. 5 (a), the type-C 
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process produces T'l^^-j {k, k; q). As shown in Fig.5(b), the type-A and type-B process 
produce T^j^^^ {k, k'; q) and Tjj^^ {k,k']q). In 7^^'^^(g), the initial spins of electron 

and hole are not flipped in the final state. On the other hand in T}}'^\k, k'; q), the 
initial spins are flipped in the final state. We note that although T]_^''^^ includes only 
one scattering channel (type-C scattering) in all the intermediate states, tJ_'^''^^ and 
q-U,lt inci^fjg different scattering channels (type-A and type-B scattering) in 
all the possible ways. As a result the number of possible diagrams corersponding to 
Fig. 6(b) increases more and more as the order of the diagram increases. 

4.2 rl^'^-^ 

We first consider T^^'^'^^k, k'; q). For example the 3rd order contribution can be 
explicitly written by 

Tll^^\k,k';q) = E Mk-ki)Mk,-k2)Mk2-k')goik,)goik,+q)goik2)go{k2+q). 

ki,k2 

(42) 

To proceed with our analysis, we need two approximations. First we neglect the 
frequency dependence of the interaction and replace J±{q) with the most dominant 
contribution J±{q) = J±{q]iuJi = 0). Then we can perform the matsubara summa- 
tions with regard to ki and k2 in (^21) and get 

rll^^ik,k';q) = E - fci)Xfc,(g;0)Jx(fci - fc2)Xfc,(g;0)Jx(fc2 - k'). (43) 

ki,k2 

Now we introduced the quantity Xkil' Q) 

XkW, Q) - -TE ^o(A: + q)g.{k + Q) = _l tanh4^-tanh -^ ^ 

We note that J±{q) has a broad maximum around the antiferromagnetic vector q* = 
(yr/g, 7r/a). Furthermore the internal wavenumber k and k' become important only 
when they satisfy the condition A; — fc' ~ {forward scattering) or fc — fc' ~ q* 
{backward scattering). The latter process is strongly driven by the antiferromagnetic 
spin fluctuations. We can see that the 3rd order process is enhanced only in the case 
of backward scattering, since — fci ~ q* and fci — ^2 ~ q* give constraint on k' as 
k' k + q*. 

This consideration gives us important results: the ladders with odd number J±{q) 
give dominant contribution to backward scattering processes, while the ladders with 
even number J±{q) give dominant contribution to forward scattering processes. This 
situations are shown in Fig. 6 Keeping this constraints in mind, we replace J±{q) 
with the averaged one Q 

J±{T) = Y.Mq,tuJi = 0) = Jo±E{7(9,^^« = 0)}^ (45) 

q q 

^To replace the non-separable momentum dependent interaction with its averaged value is in- 
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which depends on the temperature through Xo{Q] 0) the Stoner factors. This pro 

Ti,T' 
_{3) 



cedure is our second approximation. Now Tj_^^Mk, k'; q) can be approximated simply 



by 



Tl[i^\k, k'; q) = - Xo(g)^ (46) 



1^/ 



where Xoil) = ^k^k^l'^^- ^^^^ stage T-matrix depends only on q. Thus from 
now on we simply denote the T-matrix as T]_^'^^{q). 

Genelarizing the aformentioned procedure, we can get the T-matrixes for forward 
and backward process as 

^Ti,Ti.„^ _ {^rXo(g) 

i-{^xo(g)} 
Tn,n(^^ {^} xojqf 

^LBWW 7V~^\ 72' 

i-{^xo(g)} 

where FW and BW denote the case of forward and backward process respectively. 
4.3 rl^'^^ and T^^'^^ 

Next we consider TjJ^^^ and Tjj^^-^ . To get explicit expressions for them, let us 
introduce the auxiliary spin- independent T-matrix T±(n) that is defined in Fig.5(c). 
In Tl(„), all the inter-plane interactions contained in the ladder are type-B interaction. 
In reality TJJ^^^ and Tlj^^-^ contain type-A and type-B interactions in all the possible 
ways. Then n-th order ladder T^jJ^^ and Tjj^-^ can simply be related to corresponding 
71 (n) through combinatrics factor as 

'^Hn) = (49) 

n,i 

-in) 



^inV = bnT^i^) (50) 



where 

[f] 

O'n = X/^*"^2i; (51) 
1=0 



2^*^^„C2i+i. (52) 

i=0 



evitable to solve the Bethe-Salpeter equation in a closed form. In the case of Coulomb gas, the 
dielectric screening removes the singularity and as a result this approximation may work well. Wolff 
discussed the spin susceptibility of electron gas in this manner p^. 

In the present problem, the spin- fluctuation mediated interaction J±{q) has short range nature 
and doesn't have any singularity in a paramagnetic phase. Therefore in the present case this ap- 
proximation can be justified. 
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Here nCm = nl/Km — n)!m!] is a binomial coefficient. 

The above relations have already been discussed in Ref.|p4||. Now, for example, 



we derive the formula (^91). In this case n-th order term of TJJ^^^ can include the 
even number spin-flipping (Type-A) scattering. Then all the other scatterings are 
the type-B scattering. We consider the case when there are 2i {i = 0, 1, ■ ■ ■) type- A 
processses. Then there are n — 2i Type-B vertices in the n-th order ladder. Since the 
type-A vertex gives the facor Jo±/2 and the type-B vertex gives the facor Jo±/4, if 
we replace all the type-A vertices simply with the Type-B vertices to get , there 
appear the factor 2^*. Furthermore, concerning the location of the type-A interactions 
in a ladder there are „C2i ways. As a result T^J^j^ is related with '71(„) through a 

factor J:?Jo'2^\C2^. 

Now we notice that n-th order prefactors satisfy the simple relations 

a„ + bn = (1 + 2)" = 3", (53) 

and 

a„-6„ = (1-2)" = (-!)". (54) 

We can get actual feeling for these relations by writing down concrete values of mul- 
tiplication factors 
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f 02 


= 5 


/ as = 


13 


J 04 


= 41 r 


1 b2 


= 4 


'1 &3 = 


14 ' 


1 


= 40 '\ 



Using the relation, (|53|) , we can obtain 

rTT.TT I 

-±{n) + ^±(n) 



^JnV+^_;tV=3"T^H. (55) 



Therefore we can get the exphcit result for T^^'^^^ + tI^'^\ In the case of forward 
process, dominant contributions arise from 

while in the case of backward process, dominant contributions arise from 

±BW ~^ ±BW — f„ -I 2- 

Next we go back to the relation ( p4D that gives the simple relation 



_(n) - ^±(n) - y-^) ^-L(n) 

Now we can notice that 

/ 1 \nn- .... 

_(n) — ^_L(n) • 



1)"T^(„) = Tli;^. (59) 
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Therefore the relation 



directly leads to the spin-rotational symmetry relation 



Thus we can confirm again that the present treatment doesn't violate the spin ro- 
taional symmetry. 

The relations (^), ( pTf ) and ( |60D enable us to get the explicit form of tJ_^'^^ and 



■I I 



±Fwil) ■ 



1 
2 

1 
2 



1 

2 

1 
2 



{|J±(T)}'xo(g) 1 


{lJ±{T)y xoiq) 


I - {lJAT)xoiq)y 


1 - {lJ±iT)xoiq)y _ 


{lJ±{T)y Xo{q) 


{\j±{T)y Xo{q) 


} - {jJ±iT)xoiq)y 


I - {jJ±iT)Xoiq)y _ 




{lJAT)yXo{q? ' 


I - {lJ±iT)xoiq)y 


I - {\JAT)xoiq)y _ 


' {lJ±{T)y Xo{q? 1 


{lJAT)yXoiq? ' 


_i-{!J±(T)xo(g)}' 


1 - {\JAT)Xoiq)y _ 



(61) 
(62) 



(63) 
(64) 



for backward process. 
4.4 rl'""'^'' and 



By using (^Oj), (5T) and ( [UT| ) ~ (|5^, we can get the explicit form of the T-matrices 
in the singlet and triplet channel in the case of the forward and backward process as 



^singlet/ \ 

^^Fw yq) ■ 
rlJwiq) 



3 {lJ±{T)y xo{q) 1 {lJ±{T)y Xo{q) 



and 



/j-singlet 



(q) 



^triplet I 



' I - {\J±{T)xo{q)y 


'^l-{lJ±{T)xo{q)y^ 


1 {lJ±iT)y xoiq) 


1 {iJ±iT)y Xoiq) 


^ I - {lJ±iT)Xo{q)y ^ 


^l-{\J±{T)xo{q)y^ 


3 {\J±{T)y Xoiqf 


1 {lJ±iT)y xoiqf 


I - {lJ±iT)Xo{q)y 


2l-{|Jx(T)xo(g)}" 


1 {lJAT)yXoiqy 


1 {'iMT)yXo{qr 



(65) 
(66) 



1 - 



Xo 



{\J±iT)xo{q)y 



(67) 
(68) 
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In Fig. 7, we show the behavior of the singlet and triplet T-matrices as a function 
of the dimensionless quantity, J±^{T)xo{q) 1^- It follows that the inter-layer exchange 
scattering processes give rise to the strong attraction for the spin singlet pair and the 
strong repulsion for the spin triplet pair. Thus we can say that just near the intra-layer 
antiferromagnetic phase boundary the inter-layer singlet pair formation is strongly 
favored while the triplet pair formation is strongly prohibited. This tendency mainly 
comes from the dramatically enhanced repulsion between the inter-layer parallel spin 
pair, Tl^'^\q) and Tl^'^\q). We can expect that the strong repulsion for inter-layer 
parallel spin pairs induces the screening effects for the intra- and inter-layer spin 
fluctuations. We investigate this screening effects in the next section. 



4.5 Averaged Inter-layer Interaction J±{T) 

It is should be noted that the Stoner factor sensitively depends on the magnitude 
of the intra-layer Hubbard interaction. Our proposal can be justified only when 
the Stoner factor becomes considerablly large. In Fig. 8. we show the numerical 
results for the averaged inter-layer interaction Jx^ for U = 3.25t. Here the half-filling 
corresponds to /i = —0.327/!:. We fixed the chemical potential to /i = — 0.54t. By 
taking up the avarage, we can also check the non-divergence of the Stoner factor over 
the whole Brillouin zone. If the bare susceptibility touches the critical condition, 
Uxoil) = 1, at a certain temperature, Tsdw? the RPA susceptibility diverges and the 
whole theoretical backgrounds break down. 

As for the magnitude of Jo_l, we set Jo± = O.OSt. This magnitude of the bare 



inter-layre Heisenberg interaction corresponds to the result suggested in Ref . . The 
maximam of the bare susceptibility, Xo('?)) reach 0.26/t. Therefore, the divergence of 
^stngiet,tripiet ^^^-yj^.^ 0.75J7'_lXo = 1 that corrcspouds to Q.lhJ±_ ~ 3.8. As we can 
see from the result of Fig. 8, the numerical vale of 0.75j7j_ can remain finite down to 
T = with the maximum value, 0.75j7j_ ~ 3.9. If we take larger magnitude of J7j_ 
becomes divergent at unreaUistically high temperature, Tsdw- 



5 Contributions of Exchange Scattering Processes 
to Diagonal and OfF-diagonal Susceptibility 

Now we consider explicitly how inter-layer exchange scattering processees affect the 
irreducible diagonal- and off-diagonal susceptibilities. Within the present scheme, 
the inter-layer exchange processes cannot be detected directly through experimental 
probe because in our model there is no inter-layer direct carrier hopping. Instead 
these processes are inserted into the irreducible polarization. In Fig. 9 we show the 
contributions from exchange scattering T-matrices to the transverse susceptibility up 
to the 2nd order. Among them (a),(b), (c) and (d) contribute to the diagonal coun- 
terpart of the transverse susceptibility, ^xjf i although (e),(f), and (g) contribute to 
the off-diagonal counterpart, ^Xi_~ ■ Iii the appendix we shall discuss the longitudi- 
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nal counterpart, Ax^^ and Axj^, and prove the spin rotational invariance relation. 
As we shall discuss in the next section, we can notice that the neutron scattering 
experiments detect all diagrams in Fig. 9, while the NMR longitudinal relaxation 
experiments detect only diagram (a),(b),(c), and (d). Furthermore all diagrams con- 
tribute to the NMR spin echo decay rate. We note that in each diagram the forward 
and backward T-matrices are inserted correspondong to the characteristic momen- 
tum transfer. Thus we get the expressions for the correction terms to the diagonal 
counterpart as 

Ax//- = Ax^/- + Ax^/^b + Ax^/- + Ax;/>, (69) 

where 

Ax;/:a = -7(g)^r^^TT;'TT(Q)£W(g,Q), 

Q 

n„ Q 



n„ 



Q 



x{£(3)(^,g)}2 



Each indices a, b, c, ... represent the corresponding diagrams in Fig. 9. Furthermore 
the correction terms to the off-diagonal counterpart are given by 

Axr = Axt- + Axl;; + Axt-, (70) 

where 

Axte = -7(9)^r^5:r;i'J^(g)r(^)(,,Q), 

Axt7 = 7(«)'T^^[r/i,'i^(g)r7ij'^(Q + g)+ri],'J^(g)r^^^^ 

n„ Q 
nm Q 

+ ^ii'ii^(Q)'7'ii'l^(Q + + 'rlb^iQ)'rl^/,^iQ + Q)]^^'\q, Q)j^^'\-q, -Q)- 

We note that in the diagram (c), (d), (f), and (g), forward and backward counterparts 
are combined to each other. The coupling functions are given by 

£(i)(qr; iu^u Q, i^m) = TY,Y.^So{k)yGo{k + q)Qo{k + Q) 
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and 



_ y go(fc) - go(fc + q) Gojk) - go(fc + Q) 

cosh'^ ^ + xfc(0; Q) + Xfc(g; 0) - Xfc(g; Q) ^^-^^^ 



C'-^Xq; iuju Q; «a„) = ^ goik)go{k + q)go{k + Q)go(fc + Q + g) 

_ T y^y^ go{k)^gnik + q)go{k + Q)-gn{k + q + Q) 

_ -Xfc(0; Q) + Xfc(0; g + Q) + Xfc(g; Q) ~ Xfc(g; g + Q) . 



j,^,, Q; ^a^) = ^ ^0 Wao(fc + Q)g{k + Q + q) 

ga{k)go{k + Q + g) - ao(fc + Q)ao(fc + Q + g) 



= rVV. 

il7™-r;fc(Q,0) 

^ -Xfc(0; g + Q) + Xfc(Q; g + Q) /-^o\ 

^ zJ7„-r;fc(Q,0) ' ^ ^ 

where Xfc(g;Q) has already been defined by (|4^ ) and 77^(9, Q) = Cfc+g ^ ^k+Q- ^^e 
above expressions, g = {q,iu!i), Q = {Q,iQm) where cu; and flm denote bosonic 
frequencies. 

In the terminology of superconducting fluctuation, (a), (b), and (e) correspond to 
Maki-Thompson diagram|^6|], although (c), (d), (f), and (g) correspond to Aslamazov- 
Larkin diagram [Q. Maki-Thompson terms represent the coupling of spin fluctu- 
ations via the exchange scattering processes with the internal momentum transfer 
Q, although Aslamazov-Larkin term represent the coupling of the fluctuations via 
the processes with the internal momentum transfer Q and Q + q. Therefore in 
the case of the antiferromagnetic spin fluctuations, at the commensurate momen- 
tum q = q* = (tt/o, vr/a), the coupling functions in the Aslamazov-Larkin terms are 
severely reduced. Moreover in our case Aslamazov-Larkin terms inevitablly include 
7^^'^^ that is only weakly enhanced, and therefore can safely be discarded. 

Based on the above consideration, we can say that the dominant contribution to 
the diagonal susceptibility arises from the diagram Fig. 9(a), i.e. 

Ax//- ~ Ax//- . (74) 

Furthermore the most dominant contribution to the off-diagonal counterpart arises 
form the diagram Fig. 9(e), i.e. 

Axl- ~ Axl" (75) 
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Since the dominant contribution of T-matrices, 'T^±FWBwiQ)y arises from the terms 
with its external frequency Qm = 0, from now on we put Qm = in Q. Moreover the 
dominant contribution of T^lp^^B^/{flrn = 0, Q) in (|7^) and (^) arises form the term 
with q Q q*. Thus 

AxJ/-iq*;tcoi) ~ -T-f{q*;iuJifTiyi^{q*;inm = 0)C^'\q*;icui; q*;tn^ = 0), (76) 
and 

Axl-iq;tiui) ~ -T^iq;tuifrl^/^{q*;tn^ = C^^\q* ■ i^u q* ; i^l^ = 0). (77) 
Then we get the exphcit expression for couphng functions as 



6^\q*-iuJuq*-in^ = Q) = 



-2^cosh-2| ^2r7^(q*;0)Afc(g*;0) 



^ ujf + r^k{q*\^? K + r7fc(q*;0)2}2 



, (78) 



.(^)(.V.;.V.™^0) - ^E ,^(;o).F'ti;^o)^ ' '''' 
r(3)(,v...v^.^o) ^ -^E^^^^^^ (80) 

In this case C^^\q*]iuJu q*]iVLm = 0) = and thus the diagram (c), (d), (f) and 
(g) give no contribution. Ahhough we can safely say that Ax^~ gives rise to the 
negative contribution, the sign of the contribution of ^X//" is ambiguous. To settle 
this problem, it is rather instructive to investigate explicitly the behavior of the 
integrand of C^'^\q*; iuji] q*; ifi^ = 0), i.e. 

C^'\q*■^urq*■^n 0) -^'^'^'"'^ , o ^fc(g^ 0)Afc(g^ ") 
{q , q , .fi. - 0) - ^ +2^^, ^ . 

As shown in Fig. 10, the integrand is positive in almost all the region in the first 
Brillouin zone. We can see that this behavior sensitively reflects the present Fermi 
surface given by (0). Thus in the present case both of the diagram (a) and (e) give 
negative contribution and suppresses the total spin fluctuations, x + ^x/ 1 + ^X-L, in 
the low energy region. We can call this mechanism the low energy dynamical screening 
effects. From this result, we can see that the coupling of the spin fluctuations on each 
layer via the enhanced inter-layer exchange scattering processes can dramatically 
screen the low energy intra-layer spin fluctuations, x(Q')- 



6 Neutron Scattering and NMR in Magnetically 
Coupled Bi-layer Cuprates 

In the present section we show how the diagonal and off-diagonal spin-spin correlation 
functions appear in the expression of the neutron scattering cross section and NMR 
relaxation rates. 
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Neutron Scattering Cross Section 



A neutron has magnetic moment fij^ = — 7/iArS, where s is the neutron spin. The 
g- factor 7 = 1.91. This magnetic moment at the position r^v induces the magnetic 
field at the position r as B(r — tat) = Vr x (^/x^ x ^3^^) . Fluctuating electronic spin 
at the position rg interacts with this magnetic field. We note that in case of bi-layer 
compounds, is confined in layer-1 or layer-2. This interaction can be expressed by 

T^el-neu = ~(e/meC) Z]reelayer-1,2 S(re,t) • B(re — r^v)- 

The initial and final state of the electron-neutron system can be labeled by | 
kiCiEi > and | kfafEf > respectively where ka denote a wave number and spin of 
neutron and E denotes a electronic state. Then the inelestic neutron scattering cross 
section can be given by|]59[ 



da f fn^\'^ k 



dVtduj \2'nh 



\ kf 

j -rJ2 Pe^P'^, \< kiCTiEi I Hei-neu \ kfCTfEf >p 5 (energy) 



* o-io-/ EiEf 



where P^- and denote probability factor and the delta function represents energy 
conservation. By using the fiuctuation-dissipation thoerem, we get 

da _^x'J/iq,^) + ^xfiq,^) 

2^{^a(3 - qaq^) J JTf-. • (81) 



dVtduj ^ 1 — exp(— cj/T) 

From this result, we can see that both of intra- and inter-layer spin-spin correlation 
equally contribute to the inelastic neutron scattering cross section. 

NMR Longitudinal Relaxation Process 

In the case of YBCO, Hamiltonian for the spin I^j of the i-th ^^Cu (m is a 
layer index) can be written hj H = — 763H ■ (1 -)-63j^ orfe^ . j^. _|_ '}{^^_^^^^ where 753 
is the g-factor and ^^j^orfe jg ^^^q orbital contribution. The interaction between the 
^^Cu nucler spin and the neighboring electronic spins can be expressed as TCei-nuc = 
Sr^eiayer-m -^mi^'^ij^'S'^(i"j)- The hypcrfiue coupling tensor Aij is proportional to the 
amplitude of the electronic Bloch function at the nuclear site in question. In case of 
YBCO, strong hybridization between the neighboring Cu sites produces the transfer- 



hyperfine coupling in the form = dijA"^ + (5jj+a-B|61|. The first term indicates 

the core polarization, while the second term indicates the transfer hyperfine coupling 
term that is assumed to be isotropic, where i and j + a denote the nearest Cu sites. 
We note that unless there is the strong hybridyzation between the adjacent layers, we 
cannot expect the inter-layer spin correlation contributes to the Ti. 
Then the relaxation rate Ti can be obtained by K 



,82) 

^1 q ^ 

where u is the nuclear Larmor frequency. The Fourier transform of the hyperfine 

coupling for YBCO is given by ^^Aq// = A/ / + 2B{cos qxa + cos qya), ^^Aq± = A± + 2B{cosqa:a + cosqya) 
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for ^^Cu sites. We should stress that only intra-layer spin fluctuations contribute to 
the NMR longitudinal relaxation. This situation is completely different from the case 
of the neutron scattering. 

Spin Echo Decay Rate 

The spin echo detects the refocused transverse spin components after appropreate 
pulses. If the nuclear spins flip during the experiment, the spin echo decays. In 
cuprates, the spin echo decay arises from mutual spin flips caused by the indirect inter- 
action between the nuclear spins via the electronic spin fluctuations. This interaction 
is the same as the well known RKKY interaction in the case of an electron gas and can 
be written as H^^t = -{^^hf Emn {Er,A.eiayer-r.,n /^/^ Kx^'^ (r, - rif'A'^fl^^} , 
where ^x^'^i'^k — ^i) is the static non-local susceptibility. Roughly speaking, the 
Gaussian decay rate mesures the mean time of nuclear spin flips due to the indirect 
interaction. 

When we approximate the decay of the NMR spin-echo envelop by a Gaussian 

[]T0| , [6^ as e ^i^F^ the contribution from the intra-layer correlation can be written 
by 



1 



oc 



2G/ // 



E A'MxT/iQ, 0)^ - E Aq//^xTM^ 0) , (83) 



while the contribution from the inter-layer correlation can be written by 

Y ^Y.<//'^xTiQ,Of- (84) 



We note that in the inter-layer counterpart it is not necessary to eliminate the contri- 
bution from the term with = ri, since the adjacend layer are spatially separated. 
Thus the observed spin echo decay rate can be written in the form 



^2G/ V^2G/// VT2G 



1 \ 



2 



j , . (85) 



We note that both of the intra- and inter-layer spin correlation contribute to the spin 
echo decay rate. 



7 Numerical Results and Discussion 

In Fig. 11 (a) is shown the energy dependence of the retarded counterpart of the com- 
mensurate RPA dynamical susceptibility XretiQ*:^)- The retarded counterpart is 
obtained by analytic countinuation by using numerical Fade approximation |]63[| . We 



fix the chemical potential to /i = —0.54. We note that we do not treat the chemical 
potential shift due to the interaction in a self consistent manner and as a result the 
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carrier density cannot be fixed. We set Hubbard interaction U = 3.25t. As temper- 
ature becomes lower, the peak around u ~ OAQt becomes sharp. We can assign this 
peak to the dynamical nesting arising from the electron-hole excitation energy spec- 
trum Sj^^q* ~ ~ 4t(cos kxa + cos kya) at the momentum transfer q* = (vr/a, n/a). 

In Fig. 11(b), we show the energy dependence of the imaginary part of the low 
energy toto/ susceptibilities, xtet{Q*y<-^) + ^X//7et(9*' '-^) + ^xl^etl^*) that can be 
detected by neutron scattering experiments. Now x denotes the intra-layer RPA spin 
fluctuation defined by (pTj). We can see the lower side of the peak is strongly sup- 
pressed as the temperature is lowered. We can assign this suppression to dynamical 
screening due to the developed inter-layer exchange scattering processes. 

Concerning the characterisitc energy scale of the apparant-gap, we get Ug ~ OAt 
that is too large to quantitatively explain the neutron and NMR experiments. This 
situation comes from our simple treatment for the intra-layer spin fluctuations. In 
reality the intra-layer carrier propagator should be renormalized by the strongly de- 
veloped intra-layer spin fictuations. Then as the temperature is decreased, the loci of 
the peak in the dynamnical spin excitation spectrum shifts to lower energy side, due 
to the temperature dependent renormalization factor. So we can say that if we take 
into account these fictuation effects, the peak moves to the lower energy side and Ug 
can also go down to smaller value. To see the mechanism of dynamical screening in 
detail, in Fig. 12, we show separately the contributions from the diagram (a) and (e) 
in the case of T = O.OOSt. 

In our scenario, both of the diagonal and off-diagonal spin susceptibility is sup- 
pressed due to the dynamical screening effects due to the strongly enhanced inter- 
layer exchange scattering processes, and as a result the intra-layer RPA dynamical 
susceptibility that is strongly enhanced due to the intra-layer antiferromagnetic spin 
fluctuations tends to be suppressed. Thus we can get the apparent gap-like structure 
in the low energy spin excitation spectrum. 

Next we consider the NMR longitudinal relaxation rates. As we have shown, in 
our scenario, the dynamical screening effects leads to the strong suppression of the 
low energy spin excitations. Then we can expect that NMR relaxation rates should 
be suppressed. We can see from Fig. 11(b) that actually the slope at u; ^ 0, namely 
lim^^o Sx(9*) ^)/^) gradually decreases as the temperature decreases. However the 
suppression of the low energy spin excitation spectrum can be observed even in the 
simple RPA results in Fig. 11(a). In this case the suppression is due to the dynamical 
nesting. This situation shows us it is difficult to discuss NMR relaxation only in terms 
of the commensurate fluctuations. 

To clarify this point we performed the momentum integration of the contribution 
from the diagram (a) with Mila-Rice form factor |6^. The results are shown in Fig. 13. 
In this case the contribution arising from the fluctuations with q ^ q* = (vr/a, vr/a) 
becomes important since the dynamical nesting effects observed within the intra- 
layer RPA is destroyed at q q*. As a result within the intra-layer RPA, 1/TiT 
monotonically increases as the temperature decreases. By taking into account the 
bi-layer effects, however, the gap like structures survive since the internal momentum 
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of the T-matrices in the diagram (a) is independent of the external momentum and 
can be enhanced regardress of the deviation of the external momentum from the 
commensurate momentum. Thus the PSGA behavior can be obtained in the wide 
region in the momentum space and 1/TiT can detect PSGA in the present scheme. 

Finally we comment on the NMR spin echo decay rates in our scheme. As was 
mentioned in the previous section, spin echo decay rates detect the strength of the 
indirect interaction. In this case the negative sign coming from the diagram (e) is 
meaningless, and therefore the development of the contribution of the diagram (e) 
tends to enhance the spin echo decay rates. As a result, this enhancement cancels 
the screening effect arising from the diagram (a). Therefore the spin echo decay rate 
is expected to be enhanced in the low temperature region where the intra-layer spin 
fluctuations are strongly developed. This situation may be a clue to understand the 
reason why 1/TiT and I/T2G show different temperature dependence in YBCOQ 

8 Concluding Remarks 

In the present paper we proposed a possible mechanism of PSGA in magnetically 
coupled bi-layer cuprates. We shall summarize a physical picture for the pseudo-spin 
gap formation in the present scheme. The basic steps follows. 

STEP 1: In a lightly doped bi-layer cuprate, an intra-layer itinerant electronic system 
lies just near the antiferromagnetc phase boundary. 

STEP 2: The itineracy of the intra-layer system dramatically enhances inter-layer 
antiferromagnetc coupling. As a result inter-layer coupling becomes temperature de- 
pendent and developes dramatically as the temperature decreaces. 

STEP 3: Just near the magnetic phase boundary the inter-layer particle-hole ladder 
becomes very important. Dramatically enhanced inter-layer antiferromagnetic cou- 
pling induces inter-layer exchange scattering processes. This process leads to strong 
attraction for inter-layer spin singlet pair and strong repulsion for inter-layer spin 
triplet pair. Namely the strong repulsion between parallel spins is strongly enhanced 
as the intra-layer system goes nearer to the magnetic phase boundary. 

STEP 4: Strong repulsion between parallel spins lead to the dynamical screening of 
the total magnetic excitations in the low energy region. Thus PSGA can be realized 
in a magnetically coupled bi-layer cuprate. 

In the present work, we have not payed our attention to the self consistent 
renormalization effects. Now we note that the inter-layer exchange scattering pro- 
cesses do not contribute to any vertex correction, since these processes do not af- 

^Now we proceed with numerical calculations to assert this point. We will present the results 
elsewhere. 
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feet the number conservation of carriers. Thus it is sufficient to consider the self- 
energy effects. When we consider the self consistency, the diagram (a) and (b) 
in Fig. 9 should be regarded as the self energy diagram. Then the bare green 
function Qoik) should be replaced by the dressed one Q{k) = [iSn ~ ■Cfc ~ 
where the self energy is given by S(A:) = ^Y.q[T}}fw{Q) + Tl\^]^{Q)]g{k + Q). By 
putting x{k) = ?R.T.{k), z{k) = 1 — ^i^M^ can write the green function in a form 
^(^) ~ ilk) is ' where = ^^^^y^- Then the irreducible loop in RPA is replaced 

tanhlir /2'i'-tanh|. /2T 

by xiQ) = -iEfc 4k)z\k+,) i^i-ik+q+^k — ■ 

dynamical screening effects in the present context corresponds to the enhancement 
of mass-renormalization factor z{k) near the antiferromagnetic phase boundary. We 
may say that in our scenario the mass-renormalization depends sensitevely on uj and 
the low energy mass becomes heavier than the high energy mass. If this mechanism 
really takes place, PSGA is obtained again. We will report the results of sophisticated 
calculations in the forthcoming paper. 
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Appendix: Spin-rotational Symmetry Relation 

Now we show the spin rotational invariance relation in the present scheme by taking 
into account all the diagrams for the contribution to the longitudinal susceptibility. 
Necessary diagrams are shown in Fig. 14. Our goal is to show the relation 

Ax;/-(g) = 2/^xT,{q), (86) 
Axr(g) = 2Axf(?), (87) 

where 

Ax^^Jg) = \ E[Ax^n(?) - Ax^r (?)]■ 

Contribution to the longitudinal diagonal susceptibility can be obtained through the 
diagrammatic algebra shown in Fig. 15(a) and we get 

Ax;;(g) - Ax;p(g) = liqf'^ E{^i^'i^(Q) + ^ii'i^(Q)K^'n?, Q) 

Q 

+7(?)'^ i:Myu,Bw{Q)'^±Fw,Bw{Q + + 
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where 

'^±FW,BwiQ)'^±FW%wiQ + Q') = '^l^FW iQ)'^LFW (Q + Q') + "^LBW iQ)'^±BW iQ + l) 

Furthermore the algebra shown in Fig. 15(b) produces the off-diagonal counterpart 

as 

AxTiq) - Axl'-^(g) = -7(g)^^E{^-Lw(Q) - Tl'BU{Q)W'\q,Q) 

Q 

'^±FW,BwiQ)'^±BW,FwiQ + 5') '^±FW,BwiQ)'^±BW,FW 

where 

'^±FW,Bw{Q)'^±Bw!FwiQ + Q') = '^LF'W 

From now on we omitt the indices FW and BW, since these indices are not essential 
for the spin rotational symmetry. By noting the spin rotaional invariance relation for 
the exchange scattering T-matrixes, (|60D, 

we can rewrite the terms in the braces in (^) as, 

rl^'^\Q)Tl^'^\Q + q) + Tl^-^\Q)Tl-^'^-^iQ + q) - TI^-^\Q)TI^'^\Q + q) 

= Tr'^\Q)T!^'^\Q + q) + Tl^'^\Q)Tl^'^\Q + q) 

-{Tl^'^\Q) - Tl^-^^mi-rr'^HQ + q)- rl^'^^Q + q)} 

= 2{Tl^'^\Q)Tl^'^'\Q + q) + TI^'^\Q)T}^'^\Q + g)}, 

and similarly in ([89|), 

rl^'^\Q)Tl^'^\Q + q)- TI^-^\Q)TI^'^\Q + q) - TI^-^\Q)TI-^'^\Q + q) 
= {Tl^'^\Q + q) + Tl^^^\Q + q)}{Tl^^^\Q) + TI^^^\Q)] 

-Tl^''^\Q)rl^'^\Q + q)- rl^^^\Q)rl^^^\Q + q) 

= 2{Tl^^^\Q)rl^^^\Q + q) + Tl^^^\Q)rl^^^\Q + q)} 
Thus we get the results 

27(g)'^ j:{'^r'^\Q)'^±'^^iQ + q) + rl-^'^-^{Q)rr'^HQ + q)}{^^'\q, Q)}', 
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and 

Q 

Q 

These expressions are exactly twice the corresponding transverse susceptibihty, ( ]69| ) 
and (1691). Thus we obtain the spin rotational symmetry relation (pT]). 
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Figure Captions 



Fig. 1: Magnetic unit cell of YBa2Cu306. Filled and open circles denote antiparalle 
spins at planar Cu^+ sites, while shaded circles denote nonmagnetic Cu^"*" ions. 

Fig. 2: Fundamental processes induced by the inter- layer magnetic interaction; 

(a) The scattering process with spin flip (type-A), 

(b) the scattering processes between parallel spins (type-B) and 

(c) the scattering processes between anti-parallel spins (type-C). 

The straight and wavy line represent respectively the green's function of an itiner- 
ant carrier and the inter-layer antiferromagnetic interaction. Here the thick and thin 
line correspond to an carrier in the layer-1 and the layer-2 respectively. 

Fig. 3: Enhancement of the inter- layer exchange interaction vertices due to the intra- 
layer RPA spin fluctuations. 

(a) : The intra-layer transverse spin fluctuations enhances the Type-A scattering ver- 
tex. 

(b) , (c): The intra-layer longitudinal spin fluctuations enhances both of the Type-B 
and Type-C scattering vertex. The straight and wavy line represent respectively the 
green's function of an itinerant carrier and the inter-layer antiferromagnetic interac- 
tion. The dotted line represents the intra-layer Hubbard interaction. 

Fig. 4: These figures show how the intra-layer spin fluctuations in different layer 
couple to each other via the inter-layer Heisenberg interaction. The intra-layer trans- 
verse spin fluctuations can couple to each other via the type-a inter-layer scattering 
channel, as is shown in (a). On the other hands the intra-layer longitudinal spin 
fluctuations can couple to each other via both of the type-b and type-c inter-layer 
scattering channel, as is shown in (b) and (c). The shaded triangle in the right hand 
side represents the vertex enhancement due to the intra-layer spin fluctuations where 
7(5) is the Stoner factor. We can see that the inter-layer interaction is enhanced 
through the double Stoner factor and strongly enhanced. 

Fig. 5: The n-th order scattering processes for three scattering channels correspond- 
ing to the different spin dependent scattering channel. 

(a) the T-matrix that represents the electron-hole exchange scattering process be- 
tween untiparallel spin, 

(b) the T-matrix that represents the electron-hole exchange scattering process be- 
tween parallel spin, and 

(c) the auxiliary T-matrix that consists of only B-type scattering channel. 

Fig. 6: Restrictions coming from the momentum tansfer. Since the triangle vertex is 
enhanced around the momentum transfer q ~ q*, the incoming momentum (q) and 
outgoing momentum(k') of an itinerant carrier are restricted. 

Case (a): When k— k' ~ q* (backward scattering) A, only the processes including odd 
number of triangles are enhanced. 

Case (b): When k— k" ~ (forward scattering), only the processes including even 
number of triangles are enhanced. 
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Fig. 7: Behavior of T-matrices in the singlet and triplet channel as a function of 
^j'J±iT)xo- Each line represents the forward processes and backward processes 
distinctively. 

Fig. 8: The temperature dependence of the averaged inter- layer interaction Ji_{T) 
for U — 3.25t and Jo_l = O.OSi for the chemical potential ji — —0.54. 

Fig. 9: The contributions to the transverse susceptibility arising from the enhanced 
inter- layer exchange scattering process. The graph (a), (b), (c), (d) contribute to 

the diagonal susceptibility and the graph (e), (f), (g) contribute to the off-diagonal 
susceptibility. The straight and broken lines represent the propagator of the carrier 
in different layers. 

Fig. 10: The behavior of the integrand in the coupling function L^^[q*] iujf, q*] iflm = 
0) in the first quadrant of the Brillouin zone at temperature T — O.OQt and T — O.Olt. 
We set / = 1 and = -0.54t. 

Fig. 11: (a) Imaginary parts of the intra- layer RPA susceptibility for various temper- 
atures T. 

(b) The retarded counterpart of the total susceptibilities Xiet{Q*l'^) fo^ various tem- 
peratures T. We fixed the chemical potential to /j, — —0.54. We set U — 3.25t and 
Jo± = O.OSt. 

Fig. 12: The retarded counterpart of the dynamical susceptibility. We show contri- 
butions from diagram (a) and (e) separately. [RPA-I-Diagram(a)] gives the diagonal 
susceptibility that contribute to NMR longitudinalrelaxation. [RPA+Diagram(a), 

(c) ] is the sum of diagonal and off-diagonal susceptibility that contribute to the neu- 
tron scattering experiments. 

Fig. 13: NMR longitudinal relaxation rate. The dashed line represents the result 
from intra- layer RPA, while the solid line represents the result when the inter-layer 
effects are taken into account. 

Fig. 14: Diagrams that are necessary to guarantee the spin rotaional invariance rela- 
tion when we include the coupling arising from enhanced inter-layer exchange scat- 
tering processes; 

(a) for the diagonal process, and 

(b) for the off-diagonal process. 
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